Metagenomic surveillance uncovers diverse and novel viral taxa in febrile patients from Nigeria.

Effective infectious disease surveillance in high-risk regions is critical for clinical care and pandemic preemption; however, few clinical diagnostics are available for the wide range of potential human pathogens. Here, we conduct unbiased metagenomic sequencing of 593 samples from febrile Nigerian patients collected in three settings: i) population-level surveillance of individuals presenting with symptoms consistent with Lassa Fever (LF); ii) real-time investigations of outbreaks with suspected infectious etiologies; and iii) undiagnosed clinically challenging cases. We identify 13 distinct viruses, including the second and third documented cases of human blood-associated dicistrovirus, and a highly divergent, unclassified dicistrovirus that we name human blood-associated dicistrovirus 2. We show that pegivirus C is a common co-infection in individuals with LF and is associated with lower Lassa viral loads and favorable outcomes. We help uncover the causes of three outbreaks as yellow fever virus, monkeypox virus, and a noninfectious cause, the latter ultimately determined to be pesticide poisoning. We demonstrate that a local, Nigerian-driven metagenomics response to complex public health scenarios generates accurate, real-time differential diagnoses, yielding insights that inform policy.

The Prevalence of Undiagnosed Salmonella enterica Serovar Typhi in Healthy School-Aged Children in Osun State, Nigeria.

Typhoid fever remains a significant public health concern due to cases of mis-/overdiagnosis. Asymptomatic carriers play a role in the transmission and persistence of typhoid fever, especially among children, where limited data exist in Nigeria and other endemic countries. We aim to elucidate the burden of typhoid fever among healthy school-aged children using the best surveillance tool(s). In a semi-urban/urban state (Osun), 120 healthy school-aged children under 15 years were enrolled. Whole blood and fecal samples were obtained from consenting children. ELISA targeting the antigen lipopolysaccharide (LPS) and anti-LPS antibodies of Salmonella Typhi, culture, polymerase chain reaction (PCR), and next-generation sequencing (NGS) were used to analyze the samples. At least one of the immunological markers was detected in 65.8% of children, with 40.8%, 37.5%, and 39% of children testing positive for IgM, IgG, and antigen, respectively. Culture, PCR, and NGS assays did not detect the presence of Salmonella Typhi in the isolates. This study demonstrates a high seroprevalence of Salmonella Typhi in these healthy children but no carriage, indicating the inability to sustain transmission. We also demonstrate that using a single technique is insufficient for typhoid fever surveillance in healthy children living in endemic areas.

Adeno-associated virus type 2 in US children with acute severe hepatitis.

As of August 2022, clusters of acute severe hepatitis of unknown aetiology in children have been reported from 35 countries, including the USA1,2. Previous studies have found human adenoviruses (HAdVs) in the blood from patients in Europe and the USA3-7, although it is unclear whether this virus is causative. Here we used PCR testing, viral enrichment-based sequencing and agnostic metagenomic sequencing to analyse samples from 16 HAdV-positive cases from 1 October 2021 to 22 May 2022, in parallel with 113 controls. In blood from 14 cases, adeno-associated virus type 2 (AAV2) sequences were detected in 93% (13 of 14), compared to 4 (3.5%) of 113 controls (P < 0.001) and to 0 of 30 patients with hepatitis of defined aetiology (P < 0.001). In controls, HAdV type 41 was detected in blood from 9 (39.1%) of the 23 patients with acute gastroenteritis (without hepatitis), including 8 of 9 patients with positive stool HAdV testing, but co-infection with AAV2 was observed in only 3 (13.0%) of these 23 patients versus 93% of cases (P < 0.001). Co-infections by Epstein-Barr virus, human herpesvirus 6 and/or enterovirus A71 were also detected in 12 (85.7%) of 14 cases, with higher herpesvirus detection in cases versus controls (P < 0.001). Our findings suggest that the severity of the disease is related to co-infections involving AAV2 and one or more helper viruses.

Emergence and spread of two SARS-CoV-2 variants of interest in Nigeria.

Identifying the dissemination patterns and impacts of a virus of economic or health importance during a pandemic is crucial, as it informs the public on policies for containment in order to reduce the spread of the virus. In this study, we integrated genomic and travel data to investigate the emergence and spread of the SARS-CoV-2 B.1.1.318 and B.1.525 (Eta) variants of interest in Nigeria and the wider Africa region. By integrating travel data and phylogeographic reconstructions, we find that these two variants that arose during the second wave in Nigeria emerged from within Africa, with the B.1.525 from Nigeria, and then spread to other parts of the world. Data from this study show how regional connectivity of Nigeria drove the spread of these variants of interest to surrounding countries and those connected by air-traffic. Our findings demonstrate the power of genomic analysis when combined with mobility and epidemiological data to identify the drivers of transmission, as bidirectional transmission within and between African nations are grossly underestimated as seen in our import risk index estimates.

Lassa Virus Genetics.

In a pattern repeated across a range of ecological niches, arenaviruses have evolved a compact four-gene genome to orchestrate a complex life cycle in a narrow range of susceptible hosts. A number of mammalian arenaviruses cross-infect humans, often causing a life-threatening viral hemorrhagic fever. Among this group of geographically bound zoonoses, Lassa virus has evolved a unique niche that leads to significant and sustained human morbidity and mortality. As a biosafety level 4 pathogen, direct study of the pathogenesis of Lassa virus is limited by the sparse availability, high operating costs, and technical restrictions of the high-level biocontainment laboratories required for safe experimentation. In this chapter, we introduce the relationship between genome structure and the life cycle of Lassa virus and outline reverse genetic approaches used to probe and describe functional elements of the Lassa virus genome. We then review the tools used to obtain viral genomic sequences used for phylogeny and molecular diagnostics, before shifting to a population perspective to assess the contributions of phylogenetic analysis in understanding the evolution and ecology of Lassa virus in West Africa. We finally consider the future outlook and clinical applications for genetic study of Lassa virus.

Metagenomic sequencing characterizes a wide diversity of viruses in field mosquito samples in Nigeria.

Mosquito vectors are a tremendous public health threat. One in six diseases worldwide is vector-borne transmitted mainly by mosquitoes. In the last couple of years, there have been active Yellow fever virus (YFV) outbreaks in many settings in Nigeria, and nationwide, entomological surveillance has been a significant effort geared towards understanding these outbreaks. In this study, we used a metagenomic sequencing approach to characterize viruses present in vector samples collected during various outbreaks of Yellow fever (YF) in Nigeria between 2017 and 2020. Mosquito samples were grouped into pools of 1 to 50 mosquitoes, each based on species, sex and location. Twenty-five pools of Aedes spp and one pool of Anopheles spp collected from nine states were sequenced and metagenomic analysis was carried out. We identified a wide diversity of viruses belonging to various families in this sample set. Seven different viruses detected included: Fako virus, Phasi Charoen-like virus, Verdadero virus, Chaq like-virus, Aedes aegypti totivirus, cell fusing agent virus and Tesano Aedes virus. Although there are no reports of these viruses being pathogenic, they are an understudied group in the same families and closely related to known pathogenic arboviruses. Our study highlights the power of next generation sequencing in identifying Insect specific viruses (ISVs), and provide insight into mosquito vectors virome in Nigeria.

Future-proofing and maximizing the utility of metadata: The PHA4GE SARS-CoV-2 contextual data specification package.

BACKGROUND: The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatics tools and resources, and advocate for greater openness, interoperability, accessibility, and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a need for a fit-for-purpose, open-source SARS-CoV-2 contextual data standard. RESULTS: As such, we have developed a SARS-CoV-2 contextual data specification package based on harmonizable, publicly available community standards. The specification can be implemented via a collection template, as well as an array of protocols and tools to support both the harmonization and submission of sequence data and contextual information to public biorepositories. CONCLUSIONS: Well-structured, rich contextual data add value, promote reuse, and enable aggregation and integration of disparate datasets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19. The package is now supported by the NCBI's BioSample database.

Molecular profiling of the artemisinin resistance Kelch 13 gene in Plasmodium falciparum from Nigeria.

Accurate assessment and monitoring of the Plasmodium falciparum Kelch 13 (pfk13) gene associated with artemisinin resistance is critical to understand the emergence and spread of drug-resistant parasites in malaria-endemic regions. In this study, we evaluated the genomic profile of the pfk13 gene associated with artemisinin resistance in P. falciparum in Nigerian children by targeted sequencing of the pfk13 gene. Genomic DNA was extracted from 332 dried blood (DBS) spot filter paper samples from three Nigerian States. The pfk13 gene was amplified by nested polymerase chain reaction (PCR), and amplicons were sequenced to detect known and novel polymorphisms across the gene. Consensus sequences of samples were mapped to the reference gene sequence obtained from the National Center for Biotechnology Information (NCBI). Out of the 13 single nucleotide polymorphisms (SNPs) detected in the pfk13 gene, five (F451L, N664I, V487E, V692G and Q661H) have not been reported in other endemic countries to the best of our knowledge. Three of these SNPs (V692G, N664I and Q661H) and a non-novel SNP, C469C, were consistent with late parasitological failure (LPF) in two States (Enugu and Plateau States). There was no validated mutation associated with artemisinin resistance in this study. However, a correlation of our study with in vivo and in vitro phenotypes is needed to establish the functional role of detected mutations as markers of artemisinin resistance in Nigeria. This baseline information will be essential in tracking and monitoring P. falciparum resistance to artemisinin in Nigeria.

HIV-1 drug resistance and genetic diversity in a cohort of people with HIV-1 in Nigeria.

OBJECTIVE: This study was designed to provide information on the genetic diversity of HIV-1 and drug resistance mutations in Nigeria, as there is limited understanding of variants circulating in the country. METHODS: We used an advanced next-generation sequencing platform, Primer ID, to: investigate the presence of high and low abundance drug resistance mutations; characterize preexisting Integrase Strand Transfer Inhibitor (INSTI) mutations in antiretroviral therapy (ART)-experienced but dolutegravir-naive individuals; detect recent HIV-1 infections and characterize subtype diversity from a cohort of people with HIV-1 (PWH). RESULTS: HIV-1 subtype analysis revealed the predominance of CRF02_AG and subtype G in our study population. At detection sensitivity of 30% abundance, drug resistance mutations (DRMs) were identified in 3% of samples. At a sensitivity level of 10%, DRMs were identified in 27.3% of samples. We did not detect any major INSTI mutation associated with dolutegravir-resistance. Only one recent infection was detected in our study population. CONCLUSION: Our study suggests that dolutegravir-containing antiretroviral regimens will be effective in Nigeria. Our study also further emphasizes the high genetic diversity of HIV-1 in Nigeria and that CRF02_AG and subtype G are the dominant circulating forms of HIV-1 in Nigeria. These two circulating forms of the virus are largely driving the epidemic in the country.

VGEA: an RNA viral assembly toolkit.

Next generation sequencing (NGS)-based studies have vastly increased our understanding of viral diversity. Viral sequence data obtained from NGS experiments are a rich source of information, these data can be used to study their epidemiology, evolution, transmission patterns, and can also inform drug and vaccine design. Viral genomes, however, represent a great challenge to bioinformatics due to their high mutation rate and forming quasispecies in the same infected host, bringing about the need to implement advanced bioinformatics tools to assemble consensus genomes well-representative of the viral population circulating in individual patients. Many tools have been developed to preprocess sequencing reads, carry-out de novo or reference-assisted assembly of viral genomes and assess the quality of the genomes obtained. Most of these tools however exist as standalone workflows and usually require huge computational resources. Here we present (Viral Genomes Easily Analyzed), a Snakemake workflow for analyzing RNA viral genomes. VGEA enables users to map sequencing reads to the human genome to remove human contaminants, split bam files into forward and reverse reads, carry out de novo assembly of forward and reverse reads to generate contigs, pre-process reads for quality and contamination, map reads to a reference tailored to the sample using corrected contigs supplemented by the user's choice of reference sequences and evaluate/compare genome assemblies. We designed a project with the aim of creating a flexible, easy-to-use and all-in-one pipeline from existing/stand-alone bioinformatics tools for viral genome analysis that can be deployed on a personal computer. VGEA was built on the Snakemake workflow management system and utilizes existing tools for each step: fastp (Chen et al., 2018) for read trimming and read-level quality control, BWA (Li & Durbin, 2009) for mapping sequencing reads to the human reference genome, SAMtools (Li et al., 2009) for extracting unmapped reads and also for splitting bam files into fastq files, IVA (Hunt et al., 2015) for de novo assembly to generate contigs, shiver (Wymant et al., 2018) to pre-process reads for quality and contamination, then map to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences, SeqKit (Shen et al., 2016) for cleaning shiver assembly for QUAST, QUAST (Gurevich et al., 2013) to evaluate/assess the quality of genome assemblies and MultiQC (Ewels et al., 2016) for aggregation of the results from fastp, BWA and QUAST. Our pipeline was successfully tested and validated with SARS-CoV-2 (n = 20), HIV-1 (n = 20) and Lassa Virus (n = 20) datasets all of which have been made publicly available. VGEA is freely available on GitHub at: https://github.com/pauloluniyi/VGEA under the GNU General Public License.

Microbial metagenomic approach uncovers the first rabbit haemorrhagic disease virus genome in Sub-Saharan Africa.

Rabbit Haemorrhagic Disease (RHD) causes high morbidity and mortality in rabbits and hares. Here, we report the first genomic characterization of lagovirus GI.2 virus in domestic rabbits from sub-Saharan Africa. We used an unbiased microbial metagenomic Next Generation Sequencing (mNGS) approach to diagnose the pathogen causing the suspected outbreak of RHD in Ibadan, Nigeria. The liver, spleen, and lung samples of five rabbits from an outbreak in 2 farms were analyzed. The mNGS revealed one full and two partial RHDV2 genomes on both farms. Phylogenetic analysis showed close clustering with RHDV2 lineages from Europe (98.6% similarity with RHDV2 in the Netherlands, and 99.1 to 100% identity with RHDV2 in Germany), suggesting potential importation. Subsequently, all the samples were confirmed by RHDV virus-specific RT-PCR targeting the VP60 gene with the expected band size of 398 bp for the five rabbits sampled. Our findings highlight the need for increased genomic surveillance of RHDV2 to track its origin, understand its diversity and to inform public health policy in Nigeria, and Sub-Saharan Africa.

Genetic diversity and population structure of Plasmodium falciparum in Nigeria: insights from microsatellite loci analysis.

BACKGROUND: Malaria remains a public health burden especially in Nigeria. To develop new malaria control and elimination strategies or refine existing ones, understanding parasite population diversity and transmission patterns is crucial. METHODS: In this study, characterization of the parasite diversity and structure of Plasmodium falciparum isolates from 633 dried blood spot samples in Nigeria was carried out using 12 microsatellite loci of P. falciparum. These microsatellite loci were amplified via semi-nested polymerase chain reaction (PCR) and fragments were analysed using population genetic tools. RESULTS: Estimates of parasite genetic diversity, such as mean number of different alleles (13.52), effective alleles (7.13), allelic richness (11.15) and expected heterozygosity (0.804), were high. Overall linkage disequilibrium was weak (0.006, P < 0.001). Parasite population structure was low (Fst: 0.008-0.105, AMOVA: 0.039). CONCLUSION: The high level of parasite genetic diversity and low population structuring in this study suggests that parasite populations circulating in Nigeria are homogenous. However, higher resolution methods, such as the 24 SNP barcode and whole genome sequencing, may capture more specific parasite genetic signatures circulating in the country. The results obtained can be used as a baseline for parasite genetic diversity and structure, aiding in the formulation of appropriate therapeutic and control strategies in Nigeria.

Polymorphisms in Plasmodium falciparum dihydropteroate synthetase and dihydrofolate reductase genes in Nigerian children with uncomplicated malaria using high-resolution melting technique.

In 2005, the Nigerian Federal Ministry of Health revised the treatment policy for uncomplicated malaria with the introduction of artemisinin-based combination therapies (ACTs). This policy change discouraged the use of Sulphadoxine-pyrimethamine (SP) as the second-line treatment of uncomplicated falciparum malaria. However, SP is used as an intermittent preventive treatment of malaria in pregnancy (IPTp) and seasonal malaria chemoprevention (SMC) in children aged 3-59 months. There have been increasing reports of SP resistance especially in the non-pregnant population in Nigeria, thus, the need to continually monitor the efficacy of SP as IPTp and SMC by estimating polymorphisms in dihydropteroate synthetase (dhps) and dihydrofolate reductase (dhfr) genes associated with SP resistance. The high resolution-melting (HRM) assay was used to investigate polymorphisms in codons 51, 59, 108 and 164 of the dhfr gene and codons 437, 540, 581 and 613 of the dhps gene. DNA was extracted from 271 dried bloodspot filter paper samples obtained from children (< 5 years old) with uncomplicated malaria. The dhfr triple mutant I51R59N108, dhps double mutant G437G581 and quadruple dhfr I51R59N108 + dhps G437 mutant haplotypes were observed in 80.8%, 13.7% and 52.8% parasites, respectively. Although the quintuple dhfr I51R59N108 + dhps G437E540 and sextuple dhfr I51R59N108 + dhps G437E540G581 mutant haplotypes linked with in-vivo and in-vitro SP resistance were not detected, constant surveillance of these haplotypes should be done in the country to detect any change in prevalence.

Polymorphisms in Plasmodium falciparum chloroquine resistance transporter (Pfcrt) and multidrug-resistant gene 1 (Pfmdr-1) in Nigerian children 10 years post-adoption of artemisinin-based combination treatments.

The emergence and spread of Plasmodium falciparum parasites resistant to artemisinin derivatives and their partners in southeastern Asia threatens malaria control and elimination efforts, and heightens the need for an alternative therapy. We have explored the distribution of P. falciparum chloroquine resistance transporter (Pfcrt) and multidrug-resistant gene 1 (Pfmdr-1) haplotypes 10 years following adoption of artemisinin-based combination therapies in a bid to investigate the possible re-emergence of Chloroquine-sensitive parasites in Nigeria, and investigated the effect of these P. falciparum haplotypes on treatment outcomes of patients treated with artemisinin-based combination therapies. A total of 271 children aged <5 years with uncomplicated falciparum malaria were included in this study. Polymorphisms on codons 72-76 of the Pfcrt gene and codon 86 and 184 of Pfmdr-1 were determined using the high resolution melting assay. Of 240 (88.6%) samples successfully genotyped with HRM for Pfcrt, wildtype C72M74N75K76 (42.9%) and mutant C72I74E75T76 (53.8%) were observed. Also, wildtype N86Y184 (62.9%) and mutant N86F184 (21.1%), Y86Y184 (6.4%), and Y86F184 (0.4%) haplotypes of Pfmdr-1 were observed. Measures of responsiveness to ACTs were similar in children infected with P. falciparum crt haplotypes (C72I74E75T76 and C72M74N75K76) and major mdr-1 haplotypes (N86Y184, N86F184 and Y86Y184). Despite a 10 year gap since the malaria treatment policy changed to ACTs, over 50% of the P. falciparum parasites investigated in this study harboured the Chloroquine-resistant C72I74E75T76 haplotype, however this did not compromise the efficacy of artemisinin-based combination therapies. Should complete artemisinin resistance emerge from or spread to Nigeria, chloroquine might not be a good alternative therapy.

Real-time Metagenomic Analysis of Undiagnosed Fever Cases Unveils a Yellow Fever Outbreak in Edo State, Nigeria.

Fifty patients with unexplained fever and poor outcomes presented at Irrua Specialist Teaching Hospital (ISTH) in Edo State, Nigeria, an area endemic for Lassa fever, between September 2018 - January 2019. After ruling out Lassa fever, plasma samples from these epidemiologically-linked cases were sent to the African Centre of Excellence for Genomics of Infectious Diseases (ACEGID), Redeemer's University, Ede, Osun State, Nigeria, where we carried out metagenomic sequencing which implicated yellow fever virus (YFV) as the etiology of this outbreak. Twenty-nine of the 50 samples were confirmed positive for YFV by reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), 14 of which resulted in genome assembly. Maximum likelihood phylogenetic analysis revealed that these YFV sequences formed a tightly clustered clade more closely related to sequences from Senegal than sequences from earlier Nigerian isolates, suggesting that the YFV clade responsible for this outbreak in Edo State does not descend directly from the Nigerian YFV outbreaks of the last century, but instead reflects a broader diversity and dynamics of YFV in West Africa. Here we demonstrate the power of metagenomic sequencing for identifying ongoing outbreaks and their etiologies and informing real-time public health responses, resulting in accurate and prompt disease management and control.

Comparison of Light Microscopy and Polymerase Chain Reaction for the Detection of Haemoparasites in Cattle in Nigeria.

PURPOSE: Haemoparasitic diseases are among the important factors that threaten cattle health and productivity especially in the sub-Saharan region. In Nigeria, their detection using sensitive molecular techniques is scanty. This study was designed to investigate and to reevaluate the repertoire of haemoparasites of cattle in Ibadan, Nigeria with a comparative evaluation of light microscopy (LM) and polymerase chain reaction (PCR) methods. METHODS: Blood samples from 100 cattle slaughtered at Ibadan abattoirs were examined using LM and PCR techniques for haemoparasite detection. The PCR reactions using three primer sets targeting the 16S rRNA genes for Hemoplasma spp. and Anaplasma/Ehrlichia spp. and 18S rRNA genes of Babesia/Theleiria spp. were done. A few randomly selected amplicons from each set were sequenced and analysed. RESULTS: A total infection rate of 34% by LM including Hemoplasma spp. (17%), Anaplasma spp. (16%), microfilaria (5%) and Trypanosoma spp. (12%) was recorded. While, 86% positivity was recorded with PCR amplification as follows: Hemoplasma spp. (64%), Babesia/Theleiria spp. (46%) and Anaplasma/Ehrlichia spp. (5%). Comparison of LM and PCR findings showed that no LM Anaplasma spp.-positive samples and 7 out of the 17 LM hemoplasma-positive cattle were confirmed by PCR. In addition, LM led to misdiagnosis of 46 Babesia/Theleiria spp.-positive samples. Amplicon sequencing and phylogenetic analysis of Babesia/Theileria spp.-positive samples revealed Theileria velifera and Theileria annulata. In the Anaplasma/Ehrlichia spp.-positive samples, only Anaplasma marginale was characterized. Mycoplasma wenyonii, "Candidatus Mycoplasma haemobos" and Pseudomonas fluorescens like were characterized among the hemoplasma-infected cattle. CONCLUSIONS: The first report of "Candidatus Mycoplasma haemobos" and Pseudomonas fluorescens like in Nigerian cattle is herewith documented. The alarming LM misdiagnosis of haemoparasites during this study confirms its limitations as it fails to identify many parasites and emphasizes the need for inclusion of molecular techniques to improve their detection. The study also shows for the first time the high prevalence of haemotropic mycoplasma in Nigerian cattle via molecular diagnostic methods, thus indicating a strong need for the investigation of their zoonotic implications.